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We describe a numerical procedure to construct a modified velocity field that does not have any 
mean flow. Using this procedure, we present two results. Firstly, we show that, in the absence 
of mean flow, spiral defect chaos collapses to a stationary pattern comprising textures of stripes 
with angular bends. The quenched patterns are characterized by mean wavenumbers that approach 
those uniquely selected by focus-type singularities, which, in the absence of mean flow, lie at the 
zig-zag instability boundary. The quenched patterns also have larger correlation lengths and are 
comprised of rolls with less curvature. Secondly, we describe how mean flow can contribute to the 
commonly observed phenomenon of rolls terminating perpendicularly into lateral walls. We show 
that, in the absence of mean flow, rolls begin to terminate into lateral walls at an oblique angle. 
This obliqueness increases with Rayleigh number. 

PACS numbers: 47.54.+r,47.27.Te,47.52.+j 
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I. INTRODUCTION 



Rayleigh-Benard convection in a horizontal layer of 
fluid heated from below provides a canonical example 
for spatially-extended systems exhibiting pattern forma- 
tion jl] and spatiotemporal chaos [|, |j. (The lat- 
ter term refers to states that are disordered in space 
and that show chaotic dynamics in time.) In this arti- 
cle, we present results from direct numerical simulations 
of Rayleigh-Benard convection to investigate the impor- 
tance of nonlocal flow modes, generally referred to as the 
mean flow @, ^, pi [§, 61 , in the formation and dynamics 
of patterns and spatiotemporal chaos. 

In a typical Rayleigh-Benard convection experiment, 
a fluid layer is confined between two horizontal plates, 
and is thermally driven far from equilibrium by main- 
taining the bottom plate at a temperature that is higher 
than that of the top plate. As the temperature differ- 
ence is increased, the fluid undergoes an instability to a 
state in which there is motion driven by the buoyancy 
forces. When the temperature difference between the 
plates is above but near this convective threshold, a pat- 
tern comprised of patches of locally parallel convection 
rolls forms with roll diameters that are close to the depth 
of the cell; see, for example, Fig. |](a). When the tem- 
perature difference is increased, the fluid undergoes other 
instabilities that may result in the pattern developing a 
simple or chaotic time-dependence. Finally, when the 
temperature difference is increased further, spatiotem- 
poral chaotic states may a ppe ar. In particular, a state 
called spiral defect chaos [[To) is observed for Rayleigh 
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number i?~3000, Prandtl number a ~ 1, and aspect ra- 
tio r~ 16. This state is a disordered collection of spirals 
that rotate in both directions and coexist with dynamical 
defects such as grain boundaries and dislocations (see, for 
example, Fig. ||(a)). 

The following facts are known about spiral defect 
chaos: The parameter values for which it appears [Tc| , 
[EH , |l2| , [T3| , pd| ; the distributions of local roll properties 
such as the wavenumber [flEfl ; statistics of spiral and de- 
fect populations (l^, ; the mechanisms for the genera- 
tion of chaos from spatial disorder fLTflj the wavenumber 
selection mechanism for spirals [|18|, |19|; and the condi- 
tions under which spiral defect chaos transitions to other 
states. Of particular interest to this article are experi- 
ments [^0], that have observed that the spirals trans- 
form into targets when the Prandtl number is increased 
from a ~ 1 to a ^ 10 and when the Rayleigh number 
i?~3500. While this observation establishes that spiral 
defect chaos occurs only at low Prandtl numbers, it does 
not allow us to conclude which of the many dynamic 
phenomena that occur at low Prandtl numbers Q is re- 
sponsible for the formation of spiral defect chaos. 

One particular phenomenon that becomes important 
at low Prandtl numbers is the presence of mean flow. 
Mean flow is the name given to the velocity field with 
a non-zero mean over the depth of the convective cell 
that is generated by the variations of the structure of the 
convection rolls such as their curvature, amplitude, and 
wavenumber, and that in turn couples through advec- 
tion to further modify the roll structure. Its magnitude 
is approximately inversely proportional to the Prandtl 
number 

It is believed that spiral defect chaos is an effect of 
mean flow which in turn is important at low Prandtl 
numbers |l(| Q . This hypothesis has been investigated 
in numerical studies of model equations of convection. 
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For example, by coupling a mean flow-like field to the 
Swift-Hohenberg equation J23|, chaotic behavior is ob- 
served ^4|. Furthermore, when the parameter in the 
model that gives the strength of the mean flow is made 
large, spatiotemporal chaotic states akin to spiral defect 
chaos is observed ^?], Pq| . However, the above 

results are tenuous because of two reasons. Firstly, the 
Swift-Hohenberg equation describes the spatiotemporal 
behavior of a field in two dimensions, whereas convection 
is a three-dimensional phenomenon. Secondly, there are 
known limitations to the Swift-Hohenberg modelling. For 
example, it is known that the long-time dynamics of the 
Swift-Hohenberg equation may not correspond to that of 
the Boussinesq equations of convection. In particular, 
Swift-Hohenberg models exhibit spiral defect chaos as a 
transient behavior, whereas in experiments, spiral defect 
chaos is known to persist for much longer times 63]. In 
addition, it is also known that the small-scale structure 
of the mean flow at the cores of the spirals, which might 
be crucial for the persistence of spiral defect chaos, is not 
perfectly captured in the Swift-Hohenberg equation [^fj . 

Our goal in this article is to show by direct numeri- 
cal simulations of Rayleigh-Benard convection that spiral 
defect chaos is indeed a consequence of the presence of 
mean flow. In the absence of mean flow, we find that spi- 
ral defect chaotic states cease to exist, and are replaced 
by states whose statistical properties differ from those of 
spiral defect chaos. In general, studies of mean flow are 
difficult to perform in experiments, primarily because it is 
difficult to measure mean flow in an experimental setup. 
This is due to several reasons, namely that the magni- 
tude of mean flow is small (typically on the order of 1% 
of the magnitude of the velocity of the convecting rolls), 
and that it exists only in distorted and not regular pat- 
terns. To the best of our knowledge, there has only been 
one experiment that has successfully imaged aspects of 
the mean flow, but only in a simple distorted pattern 
p3| , pi| . It is not clear if such imaging techniques can be 
applied to more general and complicated patterns. Thus, 
direct numerical simulations are particularly valuable for 
the study of mean flow. 

We achieve this goal by numerically constructing a 
gedanken fluid whose velocity field is modified to have 
zero mean flow. By investigating the states that arise 
from the dynamics of this fluid and by comparing them 
with spiral defect chaos, we can infer directly the role of 
mean flow in the formation and dynamics of spiral defect 
chaos. 

Once we have the capability to remove mean flow from 
the fluid dynamics, we can apply it to the study of other 
problems. One such problem that we have investigated 
is the relation between mean flow and lateral boundaries. 
Using our numerical simulations, we have shown for the 
first time how mean flow can contribute to the com- 
monly observed phenomenon of convection rolls termi- 
nating perpendicularly into lateral walls, an observation 
that is still without much theoretical understanding. 

The remaining of this article is organized as follows: In 



Section EL we define the equations governing Rayleigh- 
Benard convection, what a mean flow is, and how it can 



be measured and eliminated numerically. In Section III 



we present results on the relation between mean flow and 
spiral defect chaos, stripe textures, and lateral bound- 
aries. In Section |VJ, we present our conclusions. 



II. DEFINITIONS 

A. Boussinesq Equations 

The evolution of a low-velocity and hence approxi- 
mately incompressible convecting fluid is governed to 
good approximation by the three-dimensional Boussinesq 
equations jl| . They are the combination of the incom- 
pressible Navier-Stokes and heat equations, with the fur- 
ther assumption that density variations are proportional 
to temperature variations and that this density variation 
appears only in the buoyancy force. Written in a dimen- 
sionless form, they are: 



_1 (d t + u • V) u(x, y, z, t) 
(dt + wV)T(x,y,z,t) 
V«u 



-Vp + V 2 u + RTzll) 
(2) 
(3) 



V 2 T, 
0. 



The field u(x, y, z, t) is the velocity field at point (x, y, z) 
at time t, while p and T are the pressure and temper- 
ature fields respectively. The variables x and y denote 
the horizontal coordinates, while the variable z denotes 
the vertical coordinate, with the unit vector z pointing in 
the direction opposite to the gravitational acceleration. 
The spatial units are measured in units of the cell depth 
d, and time is measured in units of the vertical diffusion 
time d /k, where n is the thermal diffusivity of the fluid. 
The parameter R is the Rayleigh number, defined to be 
the dimensionless temperature difference AT across the 
top and bottom plates, 



I? 



igd? 



AT, 



(4) 



where a is the thermal expansion coefficient, k the ther- 
mal diffusivity, and v the viscous diffusivity (kinematic 
viscosity) of the fluid. In this article, we will also fre- 
quently use the reduced Rayleigh number, 



R ~ Rc 
R c 



(5) 



where R c « 1708 is the critical Rayleigh number at the 
onset of convection in an infinite domain Q| . The param- 
eter a is the Prandtl number, defined to be the ratio of 
the fluid's thermal to viscous diffusivities, 



v 

a = — . 

K 



(6) 



The material walls are no-slip so that the velocity field 
satisfies 



u = 0, on all material walls. 



(7) 
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The temperature field is constant on the top and bottom 
plates: 



T 



1 \ 1 

x,y,z = =F-,i j = ±~, 



(8) 



and we assume that the lateral walls are perfectly insu- 
lating, so that 



n • VT = 0, on lateral walls, 



(9) 



where n is the unit vector perpendicular to the lateral 
walls at a given point. The pressure field p has no as- 
sociated boundary condition since it does not satisfy a 
dynamical equation. 

The influence of the lateral walls on the dynamics is 
determined by the dimensionless aspect ratio T, defined 
to be the half-width-to-depth ratio of the cell if it is rect- 
angular and the radius-to-depth ratio if it is cylindrical. 



B. Direct Numerical Simulations 

We used two different numerical schemes to solve the 
Boussinesq equations. The first is a serial second-order- 
accurate finite-difference scheme that is based on a cubic 
colocated mesh. It is highly efficient for simulating a 
rectangular cell of moderate aspect ratio. The second is 
a parallel spectral element scheme that is second-order- 
accurate in time and is able to treat more complex ge- 
ometries with arbitrary lateral boundaries. Both schemes 
were used to obtain the results presented in this arti- 
cle and were found to give good agreement with each 
other. Details of both of these schemes are available 
elsewhere |35], |3(J. For applications of these schemes 
to related problems in Rayleigh-Benard convection, see 
Refs. M, M M. 



C. Mean Flow 

When the convection pattern is made up of rolls that 
are not straight and parallel, a mean flow, slowly varying 
in the horizontal coordinates, will be set up. The im- 
portance of mean flow is that it is a nonlocal flow mode, 
and as such, affects the global behavior of the convection 
pattern even though its magnitude is small. A detailed 
derivation of mean flow can be found in Ref . . Heuris- 
tically, it can be understood as follows. When there 
are inhomogeneities in the amplitude A(x,y) and wave 
vector k.(x.y) (or equivalently, the phase <j)(x,y) where 
Vj_0 = k) of the convection rolls, a Reynolds stress will 
be generated locally from the gradients of k and A. This 
results in a flow slowly varying in the plane. In addition, 
these inhomogeneities will also induce a varying compo- 
nent p s (x, y, t) in the pressure field that is constant across 
the depth of the cell and slowly varying in the plane. The 
gradient \7±p s will then drive a global flow that, together 



with the Reynolds-stress-induced flow, distorts the con- 
vection rolls further. If we call the slowly varying flow 
uo, then we can write {221 



ad zz u D = Vj_p s 



1 

2^ 



2tt 



>u» Vuj_, 



(10) 



where the integral over the phase variable <p serves to av- 
erage out the fast modes of the integrand. The as yet 
unknown field p s can be determined via the incompress- 
ibility condition, Eq. <^j, which requires that 



,1/2 

Vj_« / dzu D (x,y,Z,t) 

J -1/2 



0. 



(11) 



Eq. (10) can then be integrated twice with respect to 
z, with boundary condition Eq. (R), to completely give 
U£>. Finally, the slow distortions, ud, advect the phase 
contours of the convection rolls, yielding an additional 
advection term in the phase equation p2[ , 



Or 



5 t + U«V_ 



(12) 



The velocity field U is called the mean flow. It is an 
average of the slow distortions over the depth of the cell, 

,1/2 

U(x,y,i)=/ dzu D (x,y,z,t)g(z), (13) 

J -1/2 

with g(z) a function that depends on the nonlinear struc- 
ture of the rolls. 

We can approximate the mean flow from our numerical 
simulations as the average over the depth of the cell of 
the slow components of the horizontal velocity, 



U(x, y,t) 



1 

2^ 



2tt 



1/2 



1/2 



dzu±(x,y, z,t). (14) 



In practice, we replace the integral over the phase vari- 
able 4> with a Gaussian filter of characteristic width 0(1) 
so that variations over short length scales are smoothed 
out. 



For the approximation of Eq. (14), and with the no-slip 
boundaries Eq. (Q), the mean flow U(x, y, t) is solenoidal: 



Vx*U = 0. 



(15) 



We will also find it convenient to use the mean flow 
stream function, Q(x, y,t), and the vertical component 
of the mean flow vorticity, lj z , defined by 



vie 



z> (V_l x U). 



(16) 



The stream function, in particular, is useful to visualize 
because it gives the stream lines and so the geometry of 
the mean flow. 
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D. Quenching Mean Flow 

We now describe a procedure to construct a modified 
velocity field that does not have any mean flow. To do 
this, we want to add to the right-hand-side of Eq. ( |l0| ) 
the negative of the source of the slow distortions, i.e., the 
depth average of the Reynolds stress, so that U£> becomes 
zero for all (x, y, z, t). In Appendix [A], we show that this 
additional term takes the form 

,1/2 

&(x,y,t) = -p I dzu-Vui (17) 

J -1/2 

with p w 1.5 a constant. We can then add $ to the fluid 
equation, so that Eq. (0) becomes 



a 1 (d t + u« V) u(x,y, z,t) = -Vp- 



V 2 u- 



RTz- 



(18) 

If 4> is introduced at time t = t q , the time needed for the 
modified velocity field u to respond to this additional 
forcing can be estimated by applying dimensional argu- 
ments on the terms in Eq. (E8h. This time scale is 0(a). 
In this paper, we consider a — 1 so we expect the mean 
flow to be quenched in a time scale of 0(1) from time t q . 

For a pattern that does not have mean flow, such as 
a pattern comprising straight parallel rolls with no de- 
fects or concentric circular rolls, the quenching proce- 
dure should leave the convective properties, such as the 
Nusselt number, of the fluid unchanged. The Nusselt 
number is the ratio of convective heat transfer to heat 
transfer that would occur by conduction alone if the fluid 
remained at rest. In Fig. [l], we show that this is true. The 
Nusselt number before (denoted by solid lines) and after 
(dashed lines) the quenching procedure, which occurred 
at t = t q = 50, are indeed the same. 



III. RESULTS 



A. Mean Flow and Spiral Defect Chaos 



Using the numerical schemes described in Section [I B , 
we evolved Eqs. (||), (||), and Jl8|) from the initial condi- 
tions 



and 



u(x, y,z,t = 0)= p(x, y,z,t = 0) = 0, (19) 



T(x,y,z,t = 0) = -z + r)(x,y,z), (20) 



where T = — z is the linear conduction profile and r\ is 
randomly chosen from a uniform distribution in the range 
[— 10~ 5 , 10~ 5 ]. We observed spiral defect chaos when the 
parameters are chosen such that the reduced Rayleigh 
number e lies in the range [0.6, 3.0], the Prandtl number 
a ss 1, and the aspect ratio lies in the range [16,30]. In 
Fig. ||(a), we show an example: a planform of the mid- 
plane temperature field T(x, y, z — 0) at time t — 500 for 
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FIG. 1: Time series of the Nusselt number N(e, t) for patterns 
comprising straight parallel rolls that have no mean flow at 
several values of e before (denoted by solid lines) and after 
(dotted lines) quenching of the mean flow, which occurred at 
time t q = 50. All data reported here are for Prandtl number 
a = 1 and a rectangular cell of aspect ratio T x — F y = 20. 
(In order to achieve straight parallel rolls, periodic lateral 
boundaries were imposed and a small sinusoidal perturbation 
in the temperature field was used as the initial condition.) 



parameters e = 1.0, a — 1, and T x = T y = 20. In general, 
the planforms we observed are qualitatively similar to 
those observed in experiments in both cylindrical |lC| ] 
and rectangular Q geometries. 



(a) 



(b) 




FIG. 2: (a) An example of spiral defect chaos observed in 
a numerical simulation using the finite-difference scheme [^BJ . 
The mid-plane temperature field is plotted at time t = 500 for 
parameters e = 1.0, a = 1, and T x = F y = 20. Dark regions 
correspond to cold sinking fluid, light regions to hot rising 
fluid. The spiral defect chaos planform is characterized by a 
disordered collection of spirals rotating in both directions and 
coexisting with dynamical defects such as grain boundaries 
and dislocations. The labels "1" to "4" are discussed in Fig. ^. 
(b) When mean flow is quenched, spiral defect chaos collapses 
to a stationary pattern of textures of stripes with angular 
bends. The planform shown here is at 10 time units after the 
quenching has been introduced to the state shown in (a) . All 
other parameters are unchanged. 

We note that the range of aspect ratios that we have 
simulated are smaller that those of past experiments 



5 



which extend up to T Z 50 0, |T|, [lj]. This limitation is 
caused by the need to integrate up to at least the hori- 
zontal diffusion time scale, t ~ T 2 , which is an estimate 
of the minimum time for thermal transients to diffuse 
over the entire cell and thus for the pattern to reach an 
asymptotic state. Because of this quadratic dependence 
on the aspect ratio, the time needed to execute the nu- 
merical schemes so that an asymptotic pattern is reached 
becomes prohibitively long for T beyond about 20 to 30 
(based on current computational resources; for details, 
seeRef. §§). 

In the rest of this Section, we report on results simu- 
lated in a rectangular cell of aspect ratio = T y = 20. 
We integrated for 500 time units, and then, at time 
t = t q = 500, invoked the forcing term given by Eq. fll7| ) 
that will quench the mean flow dynamics. In Fig. ||(b), 
we plot the mid-plane temperature field at time t = 510 
which is ten time units after the quenching of the mean 
flow has begun. (Recall that the quenching takes place 
in a time of 0(c) so the quenched state at ten time units 
should have already been asymptotic for our a = 1 state 
here.) We see that the rolls have "straightened out" 
in that they have lost their curvature and have devel- 
oped angular bends. More strikingly, the straightened 
roll patches become stationary, leaving the only dynam- 
ics in the pattern to come from the motion of defects such 
as dislocations and grain boundaries. To illustrate this 
transition from a dynamical state to a seemingly "frozen" 
one, we plot in Fig. || the time series of the rate of change 
of the temperature field at several locations in the cell. 
We see that, for t <t q = 500, the derivative dT(x,y)/dt 
fluctuates and is significantly different than zero at all 
t < t q . However, after the quenching of the mean flow 
is initiated at t = t q = 500, the derivative dT(x,y)/dt 
relaxes to approach zero in a time scale of 0(1), suggest- 
ing that all dynamics is becoming "frozen" and that a 
stationary pattern is being approached. 

We have also repeated the quenching of the mean 
flow at other Rayleigh numbers ranging from e = 
0.6 to 3.0, and for different instances of the initial condi- 
tion Eq. ( pOj ) . In all cases, we observed similar stationary 
planforms as shown in Fig. |^(b). In addition, this spiral- 
to-angular transition can be observed in the reverse di- 
rection. When the mean flow quenching is turned off at 
a later time t — 550 so that mean flow is again restored 
to the system, the angular bends develop into spirals and 
the stationary planform becomes dynamical again. Spiral 
defect chaos is fully restored [[Il|. Furthermore, the sta- 
tionary textures of stripes with angular bends can also be 
observed when the quenching is initiated at other times. 
For example, instead of initiating the mean flow quench- 
ing procedure at a time when a spiral defect chaotic state 
is already asymptotic, we have also initiated the quench- 
ing procedure immediately at the start of the simula- 



tion, t = t q = 0, again using Eqs. (|19[)- (POD as initial 
conditions. In Fig. ^(a), we show the planform after 
100 time units for the parameters e = 1.0, <r = 1, and 
= = 20. We see that it comprises patches of 
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FIG. 3: The rate of change of the temperature field dTi/dt 
versus time t for the four locations in the cell indicated in 
Fig. |^(a). Prior to quenching of the mean flow which takes 
place at time t = t q = 500, the derivative dT/dt fluctuates 
and differs from zero. After quenching, it approaches zero in a 
time scale of 0(1), suggesting that the pattern is approaching 
stationarity. 



locally straight rolls ending into each other in angular 
bends. There are no spirals present. When mean flow 
is restored at time t = 100, wc find that, after a time 
of 0(1), spiral defect chaos appears, as can be seen in 
Fig. [|(b) which shows the planform at 500 time units 
after the mean flow has been restored. 



(a) 



(b) 





FIG. 4: (a) Stationary patches of stripes with angular bends 
at time t = 100 when mean flow quenching is introduced 
at time t = 0. The parameters are e = 1.0, a — 1, and 
T x — r y — 20. (b) When the quenching is turned off at time 
t — 100 so that mean flow is restored, spiral defect chaos is 
observed. The planform shown here is at 500 time units after 
the restoration of mean flow. 



Thus, we have shown that spiral defect chaos does not 
exist without the presence of mean flow. 

Before we conclude this section, we qualitatively com- 
pare the differences between the states observed when 
mean flow is quenched and at high Prandtl numbers, for 
which mean flow is weak. (Recall that the magnitude of 



6 



mean flow is inversely proportional to the Prandtl num- 
ber.) Starting from the state shown in Fig. §(a), we 
instantaneously increased the Prandtl number for that 
state from a = 1 to cr = 10 at time t — 500. Although 
increasing the Prandtl number changes the convective 
properties of the fluid and hence the dynamics of the 
state, we nevertheless observed (see Fig. ||(a)) stripes 
with angular bends that are similar to those observed 
when the mean flow is quenched. Thus, the states ob- 
served when the mean flow is quenched and unqucnchcd 
states observed at high Prandtl numbers are similar. In 
addition, we also show in Fig. ||(b) the state observed 
when we invoke the mean flow quenching procedure after 
increasing the Prandtl number to a = 10. We see that it 
is again similar to the pattern at a = 10, suggesting that 
even at a = 10, the residual mean flow components are 
negligible. 

Finally, we note that, contrary to the results of As- 
senheimer and Steinberg [^o[ |l| , we do not observe the 
transition from spirals to targets as the Prandtl number is 
increased to a — 10. Several explanations are plausible: 
Firstly, in the Assenheimer and Stein berg experiments, 



non-Boussinesq effects are significant [20, |2l[| , whereas 
our direct numerical simulations are only for Boussinesq 
fluids. Secondly, our smaller aspect ratios may not sup- 
port the formation of targets, and that we would indeed 
see the spiral to target transition in larger aspect ratios. 
Thirdly, the transition to targets may be strongly depen- 
dent on the history of the system, in particular on the 
path (in system space) that the parameters traverse. 



(a) 



(b) 





B. Mean Flow and Pre-chaotic Stripe Textures 

At lower Rayleigh numbers near the convective thresh- 
old, the planforms observed take the form of stripe tex- 
tures rather than exhibiting spiral defect chaos. They 
comprise patches of locally parallel rolls and arcs such 
that each patch terminates at the boundaries of another 
at a different orientation, and the boundaries between 
the patches are usually populated by defects. In general, 
the stripe textures are stationary after transients, except 
for the motion of defects at the grain boundaries. In 
Fig. ||(a), we show a planform of the mid-plane temper- 
ature field at time t = 500 at e = 0.15 and a = 1 in a 
rectangular cell of aspect ratio T x = T y = 20 . 



(a) 



(b) 





FIG. 6: (a) Mid-plane temperature field at time t = 500 
for parameters e = 0.15, a = 1 and T x = T y = 20. The 
stripe texture comprises patches of locally parallel rolls and 
arcs that are stationary, (b) The state of (a) observed at 
10 time units after the mean flow is quenched. The curved 
rolls have transitioned into stripes with angular bends that 
are stationary. 

When the mean flow is quenched at time t = t q = 500, 
we observe that the stationary stripe textures remain sta- 
tionary, and that those rolls that are curved are straight- 
ened out. The resulting pattern, shown in Fig. Hb) 
which is at 10 time units after the quenching, comprises 
patches of angular structures that replaced patches of 
curved arcs. 



Nusselt Numbers 



FIG. 5: (a) The pattern observed when the Prandtl number 
is instantaneously increased from a — 1 to a = 10 comprises 
stripes with angular bends that are similar to the quenched 
patterns in Fig. 0(a). The pattern shown here is at 100 time 
units after the Prandtl number has been instantaneously in- 
creased. The parameters correspond to those of the state in 
Fig. § (b) When the mean flow is quenched for the a = 10 
state of (a), we see that the resulting pattern is qualitatively 
unchanged. Shown here is the state at 100 time units after 
the mean flow has been quenched. 



One way to quantify the changes introduced by the 
quenching procedure to a pattern is to look at its global 
convective properties, such as the Nusselt number. For 
a pattern with mean flow, its Nusselt number will be 
different than those of the unmodified velocity field of 
Eq. (|l8| ) because the latter is not a solution to the Boussi- 
nesq equations. An alternate way of saying this is that 
Eq. (|l8|), together with Eqs. (||) and (||), can be inter- 
preted as the driven Boussinesq equations with a driving 
force cr _1 $ that is turned on at time t q . Owing to this 
driving, we expect the convective properties of the fluid 
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to be stronger at time t > t q than at time t < t q . This is 
illustrated in Fig. 0. The fractional change in the Nusselt 
number AN/N caused by the introduction of the quench- 
ing of the mean flow increases with the reduced Rayleigh 
number. A best linear fit to the data yields the relation 



AN/N = (0.052 ±0.005)e. 



(21) 



Thus, for example, when e ~ 1, modifying the velocity 
field to quench the mean flow introduces a change of ap- 
proximately 5% to the averaged convective properties of 
the fluid. 
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FIG. 7: Time series of the Nusselt number N(e, i) for stripe 
textures and spiral defect chaos at several values of e before 
(denoted by solid lines) and after (dotted lines) quenching 
of the mean flow which occurred at time t q = 500. All data 
reported here are for Prandtl number a = 1 and a rectangular 



cell of aspect ratio F x 



20. 



wavenumbers of convecting spirals are "frustrated", i.e., 
they lie between two competing selection mechanisms, 
selection by focus-type singularities |4|] and selection by 
dislocations Q . These two sets of selected wavenum- 
bers, at a = 1, are denoted in Fig. || by the dashed 
and the dotted lines, respectively. We see that our di- 
rect numerical simulations produced wavenumbers (de- 
noted by the circle symbols) that lie within these two 
sets of selected wavenumbers. For comparison purposes, 
we have also included the mean wavenumbers calculated 
in a previous experiment jl^, [l4|] performed in a cylin- 
drical cell with r = 78 and a = 0.95 (diamond sym- 
bols). We see that, at lower Rayleigh numbers, the mean 
wavenumbers from our simulations agree with the exper- 
imental findings. However, at higher Rayleigh numbers, 
the wavenumbers from our simulations are smaller than 
those of the experiments. Presumably, the smaller aspect 
ratios used in our simulations mean that our wavenum- 
bers are affected by finite size effects. 
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D. Wavenumber Distributions 

In this section, we quantify the differences between the 
patterns observed with mean flow and with mean flow 
quenched by studying the wavenumber distributions. We 
compute the probability density function of wavenum- 
bers, P(k), from a time-average of the patterns. We 
used the local method discussed in Ref. [|l5| to calcu- 
late the wavenumber distributions. We have found that, 
for smaller aspect ratios T ~ 20, this method produces 
better statistics than global Fourier transform methods 
that were used in previous experiments | fio[ [l2"[ [l4|| . 
The mean of the wavenumber distribution then gives the 
mean wavenumber (k)(e) as a function of the reduced 
Rayleigh number e. 

Before we highlight the differences, we point our that 
the mean wavenumbers obtained from our numerical sim- 
ulations of spiral defect chaos lie within the Busse stabil- 
ity balloon [Q, In addition, they are also consistent 
with existing theory for the selection of wavenumbers 
in spiral defect chaos [|l8| [lj|, which suggests that the 



FIG. 8: Mean wavenumbers (k) for various reduced Rayleigh 
numbers e. The circle symbols denote wavenumbers estimated 
for spiral defect chaotic states at a — 1 and T x = F y = 20 
from our direct numerical simulations averaged over differ- 
ent random initial conditions. The cross symbols denote 
wavenumbers for states observed at 10 time units after mean 
flow is quenched. For comparison purposes, the diamond 
symbols denote wavenumbers obtained in the experiment of 
Ref. in a cylindrical cell of T = 78, a = 0.95. The dashed 
line denotes the unique wavenumber kf possessed by focus- 
selected convection at a = 1, and is represented by Eq. (p3[). 
The dotted line denotes the unique wavenumber kd selected 
by dislocations. 

For the range 0.6 < e < 1.2, the mean wavenum- 
bers of the stripes with angular bends when mean flow is 
quenched (denoted by the cross symbols in Fig. ||) appear 
to fall onto a straight line whose mathematical form can 
be obtained from a linear fit, 

(k) = (3.14 ± 0.05) - (0.16 ± 0.06) e. (22) 

This relation is consistent with the wavenumbers selected 
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by focus-type singularities at Prandtl number a — 1 [E3| , 

k f =3.117-0.13 e. (23) 

The local pattern in focus-selected convection includes 
rolls that form closed contours about a point within the 
cell. In our rectangular geometries, the four corners act 
as focus centers, as can be seen by the presence of ap- 
proximately axisymmetric roll patches emanating from 
the corners, see Fig. ||(a). In the absence of mean flow, 
the wavenumber selected therefore appears to be domi- 
nated by that selected by the focus centers (i.e., the cor- 
ners) to give a mean wavenumber consistent with that 
selected in focus- type singularities. 

Furthermore, in the absence of mean flow, the 
wavenumbers kf lie at the boundary of the zig-zag in- 
stability p^ ]. The patterns observed with mean flow 
quenched are thus dominated by lateral "zig and zag" 
bondings, leading to the stripes with angular bends ob- 
served in Figs. ||(b) and||(b). 

We have also computed the correlation length £(e) of 
the patterns as a function of e. The correlation length, 
defined here as the inverse of the standard deviation of 
the probability density function P(k), is a measure of the 
average length scale of correlated regions in the pattern. 
In Fig. ^, we show £(e) calculated for both unquenched 
patterns (denoted by the circle symbols in Fig. ^J) and 
for patterns observed when the mean flow is quenched 
(cross symbols). For comparison purposes, we have also 
included the correlation lengths calculated from a pre- 
vious experiment Q performed in a cylindrical cell 
with r = 78 and a — 0.95 (diamond symbols). We see 
that the correlation lengths for the states when mean flow 
is quenched are, on the average, about twice as large as 
those for spiral defect chaos at all values of e. In ad- 
dition, the correlation lengths for the unquenched pat- 
terns can be fitted with the power law £ oc e -1 ^ 2 , as has 
been suggested by past experiments [|l2| |l4j , and which is 
predicted by dimensional arguments to be valid at least 
near threshold. However, the same cannot be said for the 
quenched states. In fact, the data suggests that while an 
exponent of —1/2 might be fitted for e~0.7, the correla- 
tion lengths appear to have saturated at £ ~ r = 14 for 
e^O.7. This suggests that finite size effects become im- 
portant, and that, in order to obtain a better estimate of 
the scaling relation for the patterns observed when mean 
flow is quenched, we would need to use a larger aspect 
ratio. Owing to the lack of data over more decades of re- 
duced Raylcigh numbers, actual fittings to the data were 
not carried out. 



E. Curvature Distributions 

Finally, we quantify how much the quenching of the 
mean flow straightens the rolls by looking at the distri- 
bution of the local curvature Xi defined at every point in 
the planform to be the magnitude of the divergence of 
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FIG. 9: Correlation length £ vs. reduced Rayleigh number 
e. The symbols are as denned in Fig. ^| The dashed line 
corresponds to the power law £ oc e -1 ' 2 . 



the unit wave vector: 

X=|V-fe|. (24) 

A value of % = corresponds to a straight roll, whereas 
a value of x = 1 corresponds to a roll with a radius of 
curvature of unity. 

We have computed the probability density function 
P(x) f° r spiral defect chaos observed at e = 1.0, a = 1, 
and T x = T y = 20, as well as for the resulting stripes 
with angular bends observed at 10 time units after mean 
flow is quenched. In Fig. |To|, we plot the two distribu- 
tions. The curvature distribution for spiral defect chaos 
(solid line) peaks at a value of x ~ 0.1, suggesting that 
the pattern is dominated by spirals whose radius of cur- 
vature is x _1 ~ 10, consistent visually with the pattern 
shown in Fig. |2|(a). We see that this peak broadens to 
become a plateau at 5^-0.1 for the quenched state 
(dashed line) , suggesting an increase in the dominance of 
straighter rolls in the pattern. 

We observed similar results for the comparison of the 
curvature distribution for the stripe textures. In Fig. 
we show the comparisons for a state at e = 0.15. We 
see that both distributions, with (solid line) and without 
(dashed line) mean flow decrease approximately mono- 
tonically and rapidly with increasing x- Both the com- 
parisons for spiral defect chaos and for stripe textures 
suggest that the consequence of quenching mean flow is 
to straighten out the rolls. 

In addition, the distribution at e = 0.15 for the 
quenched case is higher for x~0.05 as well as for 
0.1 ~ x~ 0.4 (see the inset of Fig. [ll]), and lower other- 
wise. This suggests that another consequence of quench- 
ing mean flow is the development of angular structures 
that have large curvatures. 
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FIG. 10: The probability density function P(x) of the curva- 
ture X- The solid line is for the spiral defect chaotic state at 
e = 1.0, a = 1, and F x — T y — 20 averaged over different ran- 
dom initial conditions and times t — 400 to 500. The dashed 
line is for the stripes with angular bends observed at 10 time 
units after mean flow has been quenched. 
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FIG. 11: The probability density function P(x) of the curva- 
ture X- The solid line is for stripe textures at e = 0.15, a = 1, 
and F x = T y = 20 averaged over different random initial con- 
ditions and times t = 400 to 500. The dashed line is for the 
stripes with angular bends observed at 10 time units after 
mean flow has been quenched. Because P(x) ~ for x ~ 0-5, 
the region 0.5 ~ x ~ 1 is not plotted. The inset is a zoom of 
the region 0.1 < x < 0.4. 



F. Mean Flow and Lateral Boundaries 

In experiments where the Rayleigh number is suffi- 
ciently high, it has been frequently observed that con- 
vection rolls terminate perpendicularly into the lateral 
walls. We show in this section that mean flow generated 
by amplitude gradients near lateral walls can be used to 



explain this phenomenon, although the applicability of 
this argument rests on a number of factors, among them 
the presence of defects which affects the ability of the 
patterns to reorient themselves. 

If we call n the outward unit vector normal to the 
lateral boundary and k the wave director of the rolls, 
then we can define the wall-roll obliqueness angle as 



= arccos I k • 



(25) 



In practice, the numerical value of at a particular loca- 
tion along the lateral boundary is obtained by averaging 
Eq. ( p5| ) over a length r = 0.5 to r = 1.5, where r is the 
perpendicular distance away from that location along the 
lateral boundary. The value = 7r/2 corresponds to rolls 
terminating perpendicularly into the walls. The common 
occurrence of this value remains a phenomenological ob- 
servation, without much theoretical understanding, al- 
though it has been found |Q that is not fixed by the 
physical boundary conditions, Eqs. (0)-((9|). 

However, when rolls do not terminate perpendicularly 
at a lateral boundary, variations in the amplitude of the 
convection rolls as it decays near the lateral boundaries 
results in the generation of a mean flow. (Recall from 
Eq. ( |Io| ) that a mean flow is generated by inhomogeneities 
in the wavenumbers and amplitudes of the convection 
rolls.) The normal component (with respect to the lat- 
eral boundary) of this mean flow is cancelled by the flow 
generated from from slow pressure gradients, resulting 
in the mean flow being parallel to the lateral boundary. 
It then tends to push the rolls back to a perpendicular 
orientation. The actual calculations are worked out in 
Appendix |^. The importance of this mean flow in ensur- 
ing that the rolls terminate perpendicularly is indicated 
by observing that, in the absence of mean flow, oblique 
rolls are more prevalent. In Fig. we plot the wall-roll 
obliqueness angle averaged over the lateral boundaries for 
patterns observed at t = 500 at various reduced Rayleigh 
number e, Prandtl number a = 1 and in a rectangular cell 
of aspect ratio T x — T y — 20, with mean flow and with 
the mean flow quenched. We see that, with mean flow, 
the rolls are close to perpendicular, sa n/2. However, 
when mean flow is quenched, the rolls are more oblique, 
0~7r/2. In fact, the difference in the mean obliqueness 
between the states with mean flow and with mean flow 
quenched, A(0), increases approximately linearly with e, 



A(0) = (0.16±0.01)e, 



(26) 



as the inset of Fig. [L2| depicts. 

When mean flow is quenched, the reorientation of the 
rolls away from = ir/2 is almost instantaneous. We 
illustrate this result in Fig. |l3| for one particular reduced 
Rayleigh number, e = 1.0. In this case, the mean flow 
quenching takes place at time t = t q = 500. We see that, 
at time t = 500, the mean wall-roll obliqueness angle 
moves away from = ir/2 in a time scale of 0(1). 

The above argument that mean flow restores the rolls 
to a perpendicular orientation may not always be appli- 
cable. For example, when we performed simulations in 
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Reduced Rayleigh Number e 



FIG. 12: The mean obliqueness angle (0) as a function of the 
reduced Rayleigh number e for states with mean flow (circle 
symbols) and with mean flow quenched (cross symbols). The 
inset shows the difference between the two sets of data, A(0), 
as a function of e. 



after the mean flow has been quenched. One possible ex- 
planation might be that, in a cylindrical cell, there are 
more defects existing near the lateral boundaries and that 
these defects then pin the rolls, preventing them from re- 
orienting away from = 7r/2 when the mean flow is 
quenched. 
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FIG. 13: The change in the mean obliqueness angle (0) as 
a function of time, averaged over different random initial 
conditions. The parameters here are e = 1.0, a = 1, and 
= T y = 20. Mean flow quenching takes place at time 
t = t q = 500, so that, for 490 < t < 500, the mean wall-roll 
obliqueness angle is for a pattern whose bulk dynamics ex- 
hibits spiral defect chaos, whereas for 500 < t < 510, the bulk 
dynamics is made up of stripes with angular bends. 



a cylindrical cell of aspect ratio T = 30, we find that, 
at e = 1.0, the mean obliqueness angle (0) still remains 
close to 7r/2 when mean flow is quenched. This can be 
seen more clearly in Fig. |lj, where we show the probabil- 
ity density P(0) of obliqueness angles along the lateral 
boundaries for states observed in a rectangular cell of as- 
pect ratio T x = T y = 20 and cylindrical cell of aspect 
ratio r = 30. We see that in a cylindrical cell with mean 
flow quenched, the peak at « ir/2 is still observed 



FIG. 14: (a) Distribution of angles that rolls terminate at a 
lateral boundary in a rectangular cell of aspect ratio T x = 
T y — 20. The solid line shows the distribution for the spiral 
defect chaotic state averaged over different initial conditions 
at e = 1.0 and a = 1. The dashed line shows the distribution 
for the state with mean flow quenched, (b) The solid line 
shows the distribution for spiral defect chaos observed in a 
cylindrical cell of aspect ratio V = 30 at e = 1.0 and a = 1. 
The dashed line shows the distribution for this state but with 
mean flow quenched. 

Another instance where the above argument does not 



apply is at low Rayleigh numbers. From Eq. (B7) in Ap- 
pendix |b], the magnitude of the mean flow (Ufoc e 1 / 2 so 
that at low Rayleigh numbers, the mean flow may not be 
strong enough to reorient the rolls perpendicularly. This 
is evident in Fig. ||(a), where, at the reduced Rayleigh 
number e = 0.15, rolls are seen to terminate with an 
acute angle at the lateral walls. In this case, the pres- 
ence of a restoring mean flow can also be visualized. The 
mean flow vorticity plot corresponding to this pattern, 
shown in Fig |i~5|(a), shows the presence of strong vortic- 
ity along the bottom half of the left wall and the top half 
of the right wall. There, the restoring mean flow vortic- 
ity takes the form of long and narrow circulating "jets" 
that are about one roll size wide and several roll sizes 
long. In Fig. |l^(b) , the vorticity is plotted as a function 
of distance away from the lateral wall along the solid and 
dashed horizontal lines shown in Fig. |l5|(a). The exis- 
tence of a positive vorticity patch close to the wall and 
a negative patch further away from the wall, which to- 
gether indicates the presence of a restoring mean flow, 
agree qualitatively with the theoretical results of Fig. [l8| 
in Appendix |^. When the Rayleigh number of the state 
in Fig. |6| is increased from e = 0.15 to e = 1.0, the mean 
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flow becomes strong enough to reorient the rolls to be- 
come perpendicular to the lateral walls, and subsequently 
disappears. 




•Distance from Wall x 



FIG. 15: (a) The mean flow vorticity u> z [x,y) corresponding 
to the stripe texture of Fig. [](a) obtained using Eq. (flil). 
Light regions correspond to positive vorticity, dark regions 
to negative vorticity. The important feature in this vorticity 
map is the presence of "jet" -like structures along the bottom 
half of the left wall and the top half of the right wall, (b) 
The vorticity uj z (x) is plotted along the solid and the dashed 
horizontal lines shown in (a). The shape of ui z (x) is to be 
compared with Fig. in Appendix |b|, where a positive and 
a negative vorticity patch sets up a restoring mean flow. 



IV. CONCLUSION 

In this article, we have described a novel procedure 
to construct a modified velocity field that does not have 
any mean flow in a convecting flow. Wc have applied this 
procedure to show that spiral defect chaos does not sur- 
vive when the mean flow is quenched. Instead, a pattern 
characterized by textures of stripes with angular bends 
appears. We have also shown that the mean wavenum- 
bers of these quenched patterns approach those selected 
by focus-type singularities, which, in the absence of mean 
flow, lie at the boundary of the zig-zag instability. 

We next presented a heuristic argument on how mean 
flow can contribute to rolls terminating into a lateral 
boundary perpendicularly. We provided data to show 
that, in the absence of mean flow, rolls begin to deviate 
from a perpendicular orientation, and this obliqueness 
increases with Rayleigh number. However, the ability of 
this mean flow to restore the rolls to a perpendicular ori- 
entation may be impeded by the presence of defects that 
do not allow the rolls to reorient themselves, and at low 
Rayleigh numbers where the restoring mean flow is weak. 
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APPENDIX A: DERIVATION OF QUENCHING 
OF MEAN FLOW 

In this appendix, we derive the functional form of the 
forcing term <fr that is to be added to the fluid equation, 
Eq. (ITJ), to make the resulting fluid dynamics have zero 
mean flow. 



As mentioned in Section [I C , the mean flow comprises 
a local component generated by the Reynolds stress 
l/(27r) J Q 27r dcjyu • V^ui_ and a global component driven 
by slow a horizontal pressure gradient that is present 
in order to guarantee the incompressibility condition, 
Eq. (||). Thus, if the Reynolds stress is subtracted from 
the dynamics at all times, then mean flow will not be 
generated. We thus suggest that 



2^ Jo 



1/2 

>p / dzu'V^uj^, (Al) 

'-1/2 



where the operator p J^_y 2 dz serves as an average over 
the depth of the cell. This $ can then be subtracted 
from the fluid equation, Eq. (|lj), resulting in Eq. (18). 

We now need to evaluate the value of the constant p. 
To do this, we rewrite the equation for the slow distor- 
tions, Eq. (|i~0|), as 



ad zz u D = X7 ± p s + -!- 
Ztt 



2ji 



d<fiw Vuj 



(A2) 



Following Ref. [|22| , the Reynolds stress term near thresh- 
old takes the form 

1 f 27r 

— / #u-Vui = I(k,z)H(x,y). (A3) 
2tt Jo 



where 



and 



R(x,y) = kV ± >(kA 2 ), (A4) 



I(k,z) = w (k,z)d i 



<90o(fc,z) du>o(k,z) 



dk 2 r 



dk 2 



d z 4> Q (k,z), 



with u>o(fc, z) and 4>(k, z) the vertical profiles of the verti- 
cal velocity and the potential of the horizontal velocities, 
respectively. For systems satisfying the rigid boundary 



12 



condition, Eq. (^) , these functions are the familiar Chan- 
drasekhar functions 47 



We can then rewrite Eq. (|A2|) as 



ad zz u D = V±p s + I(k, z)R(x, y) 

fl/2 

-/' 

1-1/2 



,1/2 

i / &I(fc,^)R(x,y). (A6) 

J-l/2 



Integrating Eq. ( ]A6p with respect to z twice, and making 
use of the boundary condition, Eq. (0) , 

ctud = p(z)V^p s + J(/c,z)R(x,y) 

,1/2 

-p{z) P dzI{k,z)K(x,y), (A7) 

7-1/2 
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the Poiseuille profile, and J(k, z) the double integral of 
I(k, z) with respect to z. Employing the incompressibil- 
ity condition, Eq. (fTD), we then arrive at the equality 



12 J^*dzJ{k,z) 



/-i/2 dzIQc, z) 



(A9) 



Evaluating these integrals numerically yields p ss 1.5 for 
the rigid boundary condition, Eq. (Q). Moreover, p is rel- 
atively independent of the wavenumber k, varying from 
p = 1.4886 at k = 2.8 to p = 1.4887 at k = k c = 3.117 
and to p — 1.4886 at k = 3.4, suggesting the validity of 
treating it as a constant. 

Finally, to numerically confirm this result, we carry out 
the quenching of the mean flow, as described in Eq. (18), 
for a range of values for p, at e = 1.0 in a rectangular 
cell of T x = T y = 20. At 10 time units after effecting the 
quenching, we then measure the maximum magnitude of 
the mean flow as a function of p. We plot our results in 
Fig. Hq, where we show the maximum mean flow magni- 
tude (normalized by the maximum mean flow magnitude 
observed without quenching) vs. p for data from three 
different Prandtl numbers. We see that, when p i=s 1.5, 
the normalized maximum mean flow magnitude is indeed 
zero. 



APPENDIX B: RESTORING MEAN FLOW 
NEAR A LATERAL BOUNDARY 

In this appendix, we show that a set of straight and 
parallel rolls that are oriented obliquely at an angle 
to a lateral boundary sets up a mean flow that tends to 
restore the rolls back to being perpendicular to the lateral 
boundary. The various quantities used here are defined 
in the sketch in Fig. [l7]. 



FIG. 16: The maximum mean flow magnitude vs. various 
trial values of p. The mean flow magnitudes are normalized 
by their values at p — 0, i.e., when there is no quenching. We 
see that, when p w 1.5, the mean flow goes to zero, for all 
three Prandtl numbers. 




FIG. 17: Straight and parallel convection rolls with wave di- 
rector k terminating at a lateral boundary with outward nor- 
mal n at an angle of obliqueness O. Note that, by our defi- 
nition, < <d < 7r/2. The perpendicular distance away from 
the lateral wall is x. 



We will make the assumption that the wave vectors of 
the rolls are constant near the lateral boundary, 



k = (— kcosQ, — fesinO), 



(Bl) 



and that the convection amplitude within a correlation 
length £ of a lateral boundary is suppressed ]lj E§] : 



A(x, y) = Aq tanh 



£ cos 



(B2) 



The quantity Aq is the amplitude in the bulk. The cor- 
relation length £ = V2e~ 1/2 £ with £ = 0.385. The 
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variable x is the perpendicular distance away from the 
lateral boundary. 

Then, from the Cross-Newell equation (22), the am- 
plitude gradients near the lateral wall will result in a 
non-zero mean flow vorticity, u>, given by 



7z«V_l x [kV_L«(kA 2 )], 



(B3) 



where 7 is a constant that is inversely proportional to the 
Prandtl number a. (If we relax the assumption that the 
wavenumbers of the rolls are constant, then the compres- 
sion and dilation of the rolls as well as inhomogeneities 
in their curvatures will also cont ribute to t he m ean flow.) 
Substituting Eqs. (plh and p|) into Eq. then gives 



where the normalized mean flow vorticity, 
f(Q,x) = tan(6)sech 2 
1-3 tanh 



£ cos 



(B4) 



(B5) 



£ cos0 



is plotted in Fig. |l8j(a) for several representative values 
of 9. We see that w is positive for x/£~l, and negative 
otherwise. The currents from this vorticity pair will then 
drive the rolls back to a perpendicular orientation. 



The mean flow generated by this vorticity can also be 
easily computed. Along the lateral wall, it is given by 



|U|= 7 |k y |V^-(kA 2 ). 



(B6) 



(The component of the mean flow normal to the lateral 
wall is cancelled by the flo w co ming from the slow pres- 
sure gradient.) Using Eqs. (Bl) and (|B^), we arrive at 



\V\=2 1 A 2 k 2 C 1 g(e,x) (B7) 
where the normalized restoring mean flow magnitude in 
the direction of the lateral wall 



0(6, x) — sin(9)sech 2 — ) tanh , 

\ 4 cos By \ 4 cos B 

_ (B8) 

is plotted in Fig. |lj(b) for several representative values 

of B. 

Finally, we plot the quantity max|U| as a function 
of B in Fig. [l^. We see that the restoring mean flow 
magnitude grows monotonically from zero at B = (cor- 
responding to sets of rolls parallel to the lateral wall) to 
attain its largest value at B — > tt/2 (corresponding to sets 
of rolls perpendicular to the wall). Our analysis actually 
breaks down for |B — 7i72l ~ e 1 / 4 because modifications at 
the next order in Eq. (B2) become important |h|. 
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FIG. 18: (a) The function f(Q, x) denned in Eq. QBJ) which is 
the normalized vertical component of the mean flow vorticity. 
It is plotted here for several values of O and vs. x/£. (b) The 
function f(<d,x) defined in Eq. (B8) which is the normalized 
mean flow magnitude. It is plotted here for several values of 
O and vs. 



FIG. 19: The maximum magnitude of the mean flow U as 
a function of the wall-roll obliqueness angle O. It increases 
monotonically from zero at G = (rolls parallel to the wall). 
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